Microscopic Work Distribution of Small System in Quantum Isothermal Process 
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For a two-level quantum mechanical system, we derive microscopically the exact expression for the fluctua- 
tion of microscopic work in a multi-step non-equilibrium process, and we rigorously prove that in an isothermal 
process, the fluctuation is vanishingly small, and the most probabilistic work just equals to the difference of the 
free energy. Our study demonstrates that the convergence of the microscopic work in the isothermal process 
is due to the nature of isothermal process rather than usual thermodynamic limit condition. Our investigation 
justifies the validity of "minimum work principle" formulation of the second law even for a small system far 
from thermodynamic limit. 

PACS numbers: 05.70.Ln, 05.40.-a 



I. INTRODUCTION 



Thermodynamics usually deals with the systems of infinite 
number of degree of freedoms, in which relative fluctuations 
of the observable, e.g., energy, particle number, are inversely 
proportional to the square root of the numbers of the particles 
of the system HI. Hence for a macroscopic system consisting 
of infinite number of particles, the fluctuations are vanishingly 
small and the ensemble average can describe thermodynamic 
phenomena completely. However, concerning small systems, 
usually the fluctuations of the microscopic values of thermo- 
dynamic observable will become appreciable, and ensemble 
average alone can not longer give a complete description ^ . 
In recent years, increasing interests are drawn to the study of 
thermodynamics of small system, and the emphases are put 
on the fluctuations of the microscopic value of the observable, 
instead of their ensemble average. Some notable progresses 
have been made, examples including the Jarzynski equality 
13] HI and the Fluctuation Theorem |5|. The former connects 
the free energy difference of two equilibrium states with en- 
semble average of microscopic work in non-equilibrium pro- 
cess while the later illustrates the probabilistic "entropy de- 
crease" of a closed system within short time, or transient "vi- 
olation" of the second law. These studies shed new light on 
the understanding of non-equilibrium thermodynamical pro- 
cesses of biological motors in cells and promise important ap- 
plications to the design of small-size machines. In all these 
studies, for small systems, though fluctuations of most ob- 
servables are appreciable, there exists an exception - the work 
done during a slowest reversible equilibrium process (we use 
isothermal processes to replace slowest reversible processes 
hereafter). It has been pointed out that the fluctuation of mi- 
croscopic work done by or on a small system during a slowest 
reversible process is vanishingly small tl, JJ. Nevertheless, 
though the fluctuations of microscopic work for small systems 
in finite-time irreversible processes has been extensively stud- 
ied [6.1, and the vanishing fluctuation of microscopic work of 
classical small systems specially concerning thermodynamic 
isothermal process has been point out, to our best knowledge, 
a rigorous proof of the above result from microscopic aspect 
is still lacking, and its quantum mechanical generalization has 
not been studied yet. 



In this paper, we will investigate this problem by simu- 
lating a quantum isothermal process with infinite number of 
infinitesimal quantum adiabatic process (QAP) and quantum 
isochoric process (QIP) |8,'9','T0l. We prove rigorously from 
microscopic aspect the above result that, for a two-level sys- 
tem, the fluctuations of the microscopic work during an quan- 
tum isothermal process |10l is vanishingly small. We empha- 
size that, different from most cases in conventional statistical 
mechanics, where fluctuations vanishes in the thermodynamic 
limit, the vanishing work fluctuations for a small system in an 
isothermal process is due to the intrinsic nature of isothermal 
process. Our study also verify the universal validity of the 
"minimum work principle" formulation of the second law: it 
holds even for a small system! 



II. THE THERMODYNAMIC PROCESS IN PARAMETER 
SPACE 



We consider a two-level quantum mechanical system with 
excited (ground) states |e) (\g)) with instantaneous eigen- 
energy Ee{t) (Eg{t)) depending on time t. This two-level 
system can be modeled as a spin- 1/2 in an external magnetic 
field. It interacts with a heat bath of inverse temperature /3, 
which can be universely modeled as a collection of many 
bosons with creation (annihilation) operators a J ( a^) flTl. 
The model Hamiltonian reads [iT2l[T3l . 

H = A{t)cr, + Y^ ujqalaq + ^(A,cr_aJ + /i.e.), (1) 



where c- = \g) (e| = {ax — iay)/2 and (j^ = (|e) (e| — 
\g) {g\)/2. Initially, let the two-level system be thermalized 
to equilibrium. Then we alter the magnetic field slowly so 
that the energy level spacing A{t) slowly changes from A^ 
to Ab- During the controlling process illustrated by the 
smooth curve AB in Fig. 1, the work is done on the sys- 
tem. In the infinitely slow process, which can be alterna- 
tively regarded as a quantum isothermal process ifTOl . the 
two-level system is in the thermal equilibrium at every in- 
stant, which is described by the diagonal reduced density 
matrix ps{t) = Pe(i) |e) (e| + [1 - Pe{t)] \g) {g\, where 
Pe{t) = exp[-A(t)]/(l + exp[-/3A(t)]) satisfies the Gibbs 




FIG. 1: Schematic illustration of a quantum isothermal process 1101 
AB. Here the horizontal axis Ps is the occupation probability in the 
excited state of the two-level system, and the vertical axis indicates 
the level spacing of the two-level system. The smooth curve AB 
represents the isothermal process, whose "equation of state" can be 
expressed as A (t) — — /3^^ ln{P^^~l). The horizontal and vertical 
lines represent QIC and QAP llOJ. We can use many small QAP 
and QIP to model the quantum isothermal process. For example, we 
use a "five-step stair" path (green) A — > C — > D ■ ■ ■ — > B to 
simulate the smooth curve AB. "One-step" path (blue) and "twenty- 
step" path (orange) are also illustrated. 



distribution. It should be pointed out that, during the isother- 
mal process, there is a heat exchange between the two-level 
system and the heat bath. 

For such an isothermal process, it is difficult to calcu- 
late the microscopic work distribution directly. According 
to Ref. |8] |9] [101 . however, this process can be simu- 
lated by a series of QAP and QIP. In QAP (QIC) processes, 
there is only work done (heat exchange). Hence, using the 
changes of eigen-energies of microscopic state at instant 
t — A,C,D , we can indirectly calculate the microscopic 
work done (heat exchange) E g) SW" = Ea{C) - Ea{A) 
(dQ = Ea{D) — Ef}{C)), for a, /3 = e, g. In the parameter 
space, these QAP and QIP series processes are represented by 
the "stair" path {A — > C — > D — > • • • — > B) in Fig. 
1. When every step of the "stair" path becomes infinitesimal, 
the "stair" path becomes equivalent to the isothermal process 
AB. In this way we simulate the quantum isothermal process 
with N equal-height steps (see Fig. 1) with the small height 
A = (As — A^)/A^ where A a and A^ are the level spacings 
at point A and point B respectively. The level spacings of the 
two-level system after the (j — l)-th QIC is 



Aj- = A^ + (j - 1)A, 



(2) 



for j ~ 1, 2, • • • , A^ + 1. The initial and final point A and B 
corresponds to j = 1 and j = A^+1 respectively. When we fix 
the initial point A, and the final point B, the jump A in every 
step decrease with the increase of the step number N, and A 



approaches zero when TV becomes infinity. Obviously, when 
N — > oo, the "stair" path approaches its asymptotic behavior 
- the isothermal path (see Fig. 1). When the system reaches 
thermal equilibriums, the occupation probabilities obeys the 
Gibbs distribution defined by 



Pi 



0-/3A, [ 



=.-/3A,l-l. 



p3 

9 



pJe/3A, 



(3) 



We remark that there are three time scales in our process: 
Ta for quantum adiabatic approximation, Tc the control time 
of changing the magnetic field, hence the level spacing, and 
Tr the relaxation of the two-level system. According to Ref. 
ifTSJI , Tr is determined by the coupling strength \q (fTh. We 
consider the case that Ta ^ Tc <C Tr for a quantum adiabatic 
process where we can define the microscopic work in every 
realization of the process. 



III. MICROSCOPIC WORK DISTRIBUTION 

Having defined the "path" in the parameter space (A — Pg) 
space, we can further introduce the microscopic work and its 
corresponding probabilities for a given "path". Actually, the 
definition of microscopic work is very similar to that in Ref. 
||4||. In the above path divided into many "steps", the first step 
A — > C — > D consists of a QAP A — > C, and a QIP 
C — > D. At the beginning (the point A of Fig. 1), the sys- 
tem is initially in a thermal equilibrium state ps{A), which 
implies that the system is either in its microscopic state \g) or 
|e) with probabilities Pg and Pi respectively. We choose the 
ground state in the energy reference point so that the micro- 
scopic energy E{A) of the system at initial pint A can take 
E^{A) = Aa or Eg{A) = 0, with probabiHty P} and 1 - P} 
respectively. In the first QAP A — > C, the system remains 
in its microscopic state \g) (|e)) if the system is initially in its 
microscopic state \g) (|e)). As there is no heat exchange in the 
QAP, the work done by external controller is just the change of 
the microscopic energy Wa — Ea{C) — Ea{A) for a — e, g. 
Correspondingly the work done during A — > C can be either 
Ac — A^ or with probabilities P^ or 1 — P^ respectively. 
This also agrees with the definition of work in quantum me- 
chanical system: work is associated with the change of the 
level spacing |9, 10, 14 1. 

After the QAP, a quantum isochroc process C — > D (see 
Fig. 1) follows. Here, there is no work done according to the 
definition of work in quantum mechanical system |9, 10, 14|, 
because there is no change in the eigenergies. Nevertheless, 
there is heat exchange between the system and the bath. The 
QIP last long enough i':^ Tr) so that the system can reach 
thermal equilibrium with the heat bath. After a thermolization 
for long-time, the two-level system reach thermal equilibrium 
with the heat bath again ([3]) at instant D indicated in Fig. 1. 
Then a second step D — > E — > F begins. Similarly, the 
microscopic work or A3 — A2 is done in this step with prob- 
abilities 1 — Pg or P^. The microscopic work done and their 
probabilities for the remaining steps can be obtained through 
a similar analysis. Because in every QIP, the system is inde- 
pendently thermalized by the heat bath, then there should be 
no correlations of the probabilities distributions in every two 



neighbor steps, or alternatively, this process is Markovian pro- 
cess. Hence, the total microscopic work done after the iV-step 
is a sum of microscopic works done in all steps and the joint 
probabilities for the iV-step as a whole is the product of that 
of all steps. 

For a special example that the microscopic work done dur- 
ing the whole process is W — NA where A is that for each 
QIP step, the joint probabilities for the system keeping in |e) 
in every QIP is P [NA] = P^P^ ■ ■ ■ P^ ■ The more gen- 
eral case with microscopic work W = {N — k) A corre- 
sponds to a microscopic process, in which k out of N QIPs 
ends with the system in its microscopic state \g). The prob- 
ability P{k) :— P \{N — k) A] with the microscopic work 
W — {N — A;) A in the iV-step path is given by the following 
eqution: 



N 






'^-l (,I3A, 



pt3{AA+lA) 



J=0 



g/3(/+l)A _ I 



(4) 



To prove the above result, we first consider the case with 
k — 1. For this case, there is one and only one out of the 
N QIPs, in which the system ends up in the microscopic state 

\g). Then the corresponding probability can be caculated as 

P(l)=(l-Pi)p2...pf + Pi(l-p2)...p^^ + ...+ 



plp2 



i-pn = [n7=iPi 



p(i) 



o/3A, 



2^X1=1 



g-/3Axi Qj. 




(5) 



That means thre Eq. (J4| holds for k — 1. Similarly 
we can check the case with k ~ 2. For this case, there 
are two out of the N QIP, in which the system ends up 
in the microscopic state \g). Hence its probability can 
be expressed as P(2) = (l - P^) (l - Pe) P! ■ ■ ■ Pf ' 
(l-Pi)P,2(l-p3)...pf + . 

(1 - pn = (uu ^) ^''^^ ^Z^i ^ir^i 



P^P^ • • • (1 - P. 



JV-l^ 



or 



N 



pm= \Y[Pi 



[e/3A, 



g/3A^l rg/3A, 



„/3(Aa+A)] 



= 1 



(e/3A _ I)(e2/3A _ 1) 



(6) 

Hence Eq. (J4|i also holds for k = 2 case. In general, for an 
arbitrary k, the corresponding probability can be expressed as 



N 



pik) = n Pe 



„fe/3Ai 



X{k), 



(7) 



\3 = ^ 



Where x{k) = E^^, EllZl ■ ■ ■ ^^Z^' e-'^^(-+-+-). 
As x{k), {k = 1,2,- ■ ■ , N) satisfy the following relation 



we can use the complete induction method to prove that the 
x{k) can be generally expressed as 



k-l g-/3JVA \pj3NA _ g/3iA-| 

X{k) ==11" 



al3{l + l)A _ I 



(9) 



Substituting Eq. (|9]l into Eq. (|7]i, we obtain Eq. Bjl. Hence, 
by now we prove the general result given by Eq. (Hi 



IV. MOST PROBABILISTIC DISTRIBUTION AND 
FLUCTUATION 

The above equation (Hll can result in the main conclusion 
in this paper. From the above microscopic work distribution 
function (|4|, we obtain the ratio R{k) = P(fc + 1)/P(fc) of 
distributions for two close microscopic work, i.e.. 



R{k) 



aPAE 



^0AA+kpA 



gl3{k+l)A „ I 



(10) 



Let k maximaze the probability distribution P(fc) for the mi- 
croscopic work [N - {k+ 1)]A . Then P(fc) > P{k ± l),or 
R{k) < 1 or R(k - 1) > 1. For very large k, R{k) ~ 1 that 



fcA = - In 



1 + cxp[/3Ab] 



exp[/3(AB-A^)/iV]+exp[/?A, 



(11) 



In the large N limit, the above equation determines the micro- 
scopic work W = {N — fc) A with most probabilistic distribu- 
tion 



1 / 1 + e'^'^^ 



(12) 



which is just the free energy difference AFab = Fb — Fa 
where P, = ln[l + exp(/3Aj)]//3 for j = A,B 

Next let us give a heuristic analysis of the dispersion of 
the work distribution (J4]i. Because all steps in the "stair" 
path are independent with each other, thus the whole process 
can be regarded as Markovian. So the variance of total mi- 
croscopic work done during the whole process equals to the 
sum of variance of local microscopic work in every step, i.e., 

(WXb) - k'^AB? = Ef=i((l^/) - (W^.)'), where W, is 
the microscopic work done during the jth QAP, and the local 
fluctuations 



(W^,)-{W,f^A\Pl^{Plf 



(13) 



for different j are similar. Here A is inversely proportional to 
N , and ^y/ ab) being independent of N , the relative variance 
of Wab is inversely proportional to 



fc-i 



x(fc) = E 



rii^^ H-'")x(fc-) 



vJ=i 



(8) 



Hi 



- e-J'^ 



g-fc(Ar+l)A _g-fcA 



'(w^: 



ab) 



{W^ 



AB) 



1 



{W^ 



ABl 



(14) 



We numerically plot the work distribution function (see 
Fig. 2) based on the above analytical result (J4|i to test the above 
analysis. Here we choose the step number A^ from 1 to 10000. 




FIG. 2: Microscopic work distribution of an TV-step "stair" pro- 
cess. The horizontal axis indicates the possible microscopic work 
ranging from to As — A a, and the vertical axis is their prob- 
abilities. Here, exp(-/3AA) = 1/2, and exp(-/3As) = 1/3. 
The steps are chosen to be Af = 1, 5, 20, 100, 1000, and 10000 
respectively. "Path" corresponding to A'^ = 1, 5, 20 are given 
in Fig. 1. From these figures it can be inferred that when A'^ is 
small the process is irreversible, and the fluctuation is apprecia- 
ble. The relative fluctuation of the microscopic work vanishes when 
N —> oo, or the fluctuation of an isothermal process approaches 
zero. Besides, the most probabilistic work from the (numerical) fig- 
ures W = 0.29(ln 3 - In 2)fcsT agrees well with the (analtical) free 
energy difference AFab = [ln(l + 1/2) - ln(l + l/3)]fcflr. 



V. MINIMUM WORK PRINCIPLE FOR A TWO-LEVEL 
SYSTEM 

As we have mentioned before, for small systems and within 
short time, the formulation "entropy never decrease for a 
closed system" of the second law may be transiently "vio- 
lated" probabilistically due to appreciable fluctuations |5|. A 
straightforward question is: will the other formulations of the 
second law, e.g., the minimum work principle Q [151 . also 
be transiently "violated" probabilistically for small systems? 
The "minimum work principle" states that "when varying the 
speed of a given process for an initially equilibrium system, 
the work is minimal for the slowest realization of the process" 
Q [TSJI . In the following we will test the validity of "mini- 
mum work principle" for a two-level system by utilizing the 
formula (Hli we derived above. The average work over all pos- 
sible realizations for a given A^-step path can be expressed as 



N I N 
fc=0 \j = l 




„/3(A^+iA) 



g/3(i+l)A _ I 




<W>_N 
0.1354- 



0.130 



0.125 



0.120 



LnN 



For iV = 1, the "stair" path becomes a "one-step" path con- 
sists of an QAP and an QIP (see Fig. 1). The microscopic 
work corresponding to the "one-step" path is either A^ — A^ 



or with the probability P{W 



A^) 



Pi 



P(W = 0) = 1 — Pg . In the above figures, we choose 
exp(— /3Ayii) = 1/2, (Pg = 1/3), and the numerical result 
agrees well with our analysis. For N ~ 5 (see Fig. 1), the 
possible microscopic work can be M^ = {{Ab — Aa)/5, i = 
0, 1, 2, • • • ,5. The numerical result indicates vanishing prob- 
ability for W = Ab - Aa- For TV = 20 (see Fig. 1), the 
numerical result show even more vanishing probabilities of 
microscopic work. That is, the dispersion (fluctuation) of mi- 
croscopic work decrease with the increase of N. Actually, 
from the above numerical figures, it is not difficult to find 
that the dispersion of the microscopic work distribution is in- 
versely proportional to the square root of N. For example, the 
dispersion for N ~ 100 is ten times that for N = 10000 case 
(see Fig. 2). Hence, numerical results agrees well with our 
heuristic analysis and both they verified our main result, when 
N — > oo, the fluctuations of microscopic work vanishes. 



FIG. 3: Averaged work (W) j^ as a function of A'^ 1 15 i. The steps 
Af chosen here are 5" = 1,5^ = 5,5^ = 25,5^ = 125,5* = 
625, 5^ = 3125, and 5^ = 15625. It can be seen that the av- 
eraged work is a monotonically deceasing function of A'. In the 
one-step path (A'^ = 1), the averaged work equals to (W)-^ = 
(In 3 - ln2)/3 ^ 0.135155^3?; In the 15625-step path, the av- 
eraged work equals to {^1^)56 = (ln3 - ln2)/3 ^ 0.11784fesr, 
which is very close to its asymptotic value AFab = [lii(l + 1/2) — 
ln(l + l/3)]fcsr Ki 0.117783fcflr. Thus, it can be inferred that the 
"minimum work principle" still holds for a two-level system. 

In Fig. 3 we plot the averaged work {W)j^ as a function of 



A^ ( 15 I. It can be seen that for the two-level system, {W)j^ 
is a monotically decreasing function of N (time t), and when 
TV — > 00, {t — > 00), the averaged work {W)j^ approaches 
an asymptotic value, and its minimum value - the difference 
of the free energy. Thus, from the numerical result it can be 
inferred that the "minimum work principle" still holds for a 
two-level system. 

The above proof of minimum work principle can be alter- 
natively understood in the following way. From the above 
analytical and numerical result, we observed that the fluctua- 
tion of microscopic work in an isothermal process vanishes. 



and then the work of the most probabiHstic distribution equals 
to the difference of the free energy W = AF. According to 
Ref. ID, (Wii-re) ^ AF, where (Wjrre) is the average work 
done during an irreversible process. Combining the two re- 
sults, we have (VFirro) ^ W. Thus we proved the minimum 
work principle for small system. 



VI. DISCUSSION AND CONCLUSION 

Before concluding this paper, we would hke to empha- 
size the following points: First, the technique of simulating 
isothermal processes with adiabatic processes and isochoric 
processes are important to our proof, which enables us to es- 
tablish the connection between large time limit and large N 
limit. Second, the calculation of exact expression of micro- 
scopic work in our paper is non-trivial because the work con- 
tributions in the different steps are not identically distributed. 
Hence, it is different from the law of large numbers, with time 
as the large number Q. Third, we proved the "minimum 
work principle" formulation of the second law stands for even 
small system, though other formulations may be transiently 
"violated" probabilistically jSl- This is not surprising because 
"minimum work principle" concerns infinite-long-time pro- 
cesses, which has no contradiction with the transient "viola- 
tion" of the second law for small systems predicted by the 
Fluctuation Theorem. Actually, the Fluctuation Theorem does 
not constitute real violation of the second law, which is a sta- 
tistical law and holds when averaged over different realization 
of the process. Fourth, the isothermal process is reversible, 
but the finite N "step path" is irreversible, due to the QIP 
(thermolization) is irreversible. We can thus expect that the 
work dissipation [|6]|2l for the finite N step path will be finite 
and will decrease with the increase of N, and finally vanishes 
when N approaches infinity. 

In summary, by simulating an quantum isothermal process 
with infinite many infinitesimal QAP and QIP, we obtain the 



analytical expressions of microscopic work distribution in an 
isothermal process. Through both analytical and numerical 
analysis, we rigorously verify that the fluctuations of micro- 
scopic work distribution vanishes even for a small system in 
an isothermal process. This result is different from the usual 
fluctuations in statistical mechanics, e.g., the energy fluctu- 
ation and particle number fluctuation in canonical ensamble 
and grand canonical ensamble, where the fluctuations of en- 
ergy and particle nubmers approaches zero when the system 
approaches thermodynamic limit (particle number approaches 
infinity Np — > oo). Here, however, even for single particle 
system, we microscopically demonstrate the vanishing of mi- 
croscopic work fluctuation. Because N — > oo is a must to 
simulate an isothermal process, we conclude that the vanish- 
ing of microscopic work fluctuations is due to the intrinsic 
nature of isothermal process, rather than the thermodynamic 
limit of the system size. We also prove that for a small system, 
the "minimum work principle" formulation of the second law 
holds though other formulations maybe transiently "violated" 
probabilistically. Finally we would like to point it out that our 
result is universal and does not depend on the specific model 
used here, because the technique of simulating the isothermal 
process with the isochoric process and the adiabatic process 
can be applied to any systems. Generalizations of our current 
discussion to other models will be given in the future. 
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